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ABSTRACT 

New public numerical code for fast calculations of the cosmological recombi- 
nation of primordial hydrogen-helium plasma is presented. The code is based 
on the three-level approximation (TLA) model of recombination and allows 
us to take into account some fine physical effects of cosmological recombi- 
nation s imultaneously with using fudge factors. The code can be found at 
http: / /www.ioffe.ru/ astro/QC/CMBR/atlant / atlant .html| 
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1 INTRODUCTION 

Cosmological recombination is one of the key pro- 
cesses in the early Universe. It determines the epoch of 
decoupling of radiation from matter and thereby deter- 
mines the epochs at which baryonic matter can start to fall 
into gravitational potential wells created by clustering cold 
dark matter (CDM). Afterwards these CDM+baryonic mat- 
ter clouds develo p into non-relativisti c gravity bound sys- 
tems l i ke galaxies (iPeeblej (Il965l.ll968l 'l ; iDoroshkevich et akl 
l|l967l '): lMa fc BertschingeJ l| 19951 ) 1. The sizes of these proto- 
objects also depend on cosmological recombination via the 
kinetics of divergence of radiation and matter temperatures 
which determines the critical Jeans lengt h. Cosmological re- 
combi nation affects primordial chemistry (|Dalgarno fc Leppl 
j 19871 )) and correspondingly rate of radiative cooling of col- 
lapsing clouds (via emission in resonant lines of molecules 
which depends on the abundances of primordial molecules) . 

From observational point of view the cosmological re- 
combination is also very important process because its ki- 
netics determines the positi o n of last s catter ing surface 
l|Sunvaev fc ZeldovichI (|l970l ). IHu et all (|l995l )). This in 
turn affects the cosmic microwave background radiation 
(CMBR) anisotropy which is one of the main sources of in- 
formation about evolution of the early Universe, its com- 
position and other properties. A great number of exper- 
iments on CMBR anisotropy have been carried out in 
the last thirty years (Relikt-1 1983, COBE 1989, QMAP- 
Toco 1996, BOOMERanG 1997-2003, MAXIMA 1998-1999, 
WMAP 2001-present day, and many others). Treatment 
of the results of these experiments demands clear under- 
standing of cosmological recombination physics. That is 
why many efforts have been made for theoretical inves- 



tigation of cosmological recombination and development 
of applied numerical codes for modelling of this process 
(e.g. Z cldovich et al. (1968), Peebles (1968), Matsu da et all 
1 19691). I Jones fc Wvsel (|l985l ). lGrachev fc DubrovichI l|l99ll ). 
ISeager et all (|l999l )). Increasing accuracy of experiments in 
the last decade leads to increase of efforts of theorists in 
the study of details of cosmological re combin a tion ( s ee e., 
Dubrovich fc Grachev! 1200 5 ). Chlu ba et all (120071, [201 



Chluba fc SunvaevI ll2009al. l2010bi). iRubino-Marti^n et all 



ll2008l).ISwitzer fc Hiratal (l2008bl ). iHirata fc Switzeij (|2008l ). 
lAli-Hai'moud et al.l ( 2010l )) and perfection of numerical 
codes. In the light of soon releases of experimental data from 
Planck mission this task becomes more and more important 
and urgent (for overview of efforts of investigators to find 
exact recombin ation scenario and to estim a te remain uncer- 
taintie s see e.g. ISunvaev fc Chlubal (|2009l ): [Shaw fc Chlubal 
(l201ll )). Today for successful treatment of Planck data the 
minimal required accuracy of numerical codes evaluating 
cosmological recombination is about 0.1% (on free electron 
fraction) for the epoch of hydrogen recombination and 1% 
for the epoch of helium recombination. Desirable accuracy 
is about 0.01% and 0.1% correspondingly. 

The next possible step of the investigations of 
the early Universe in the epochs z — 800 — 10^ 
is the experimental study of CMBR spectral distor- 
tions originated from cosmological recombination of hy- 
drogen and helium. Such experiments would be pow- 
erful sources of information about history of the Uni- 
verse in these epo chs. Thanks to the numerous theoretical 
works i n this field jDubrovichl (Il975l). Lvubarskv fc SunvaevI 



( 19831) :lFahr fc LochI (ll99lh:lRvbicki fc dell Anton io! (^1993 
Boschan fc Biltzingeij ll 19981 ) : lDubrovich fc Grac hev ( 2004|; 
Rubino-Martm et al. (|2008l ) and references therein) one can 
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understand parameters of these spectral distortions clearly 
enough. Of course the cosmological recombination is an in- 
tegral part of modelling of these di stortions. Note th a t sorti e 
simple metho ds for this problem (jBernshtein et all |l973); 
iBurginI l|2003h ') demand the knowledge of derivatives of ion- 
ization fractions of hydrogen and helium, so calculated ion- 
ization fractions should be smooth as possible (i.e. without 
sharp numerical features) . In spite of the fact that detection 
of CMBR spectral distortions by cosmological recombina- 
tion is impossible today, the rapid progress of measurement 
equipment for CMBR anisotropy allows us to hope that such 
experiments will be possible in the near future. 

Thus the main aim of this work is to present nu- 
merical code covering investigations of cosmological recom- 
bination widely as possible within simple model which 
we used. Our code called atlant (advanced three level 
approximation for numerical treatment of cosmological re- 
combination) may be useful for regular calculations of free 
electron fraction for treatment of CMBR anisotropy data, 
further investigations of cosmological recombination, theo- 
retic al pr edictions of new observational cosmological effects 
fe.g. lDubr ovich ct al. (2009); Grachev & Dubrovich (2010)), 
an d comparison with other numerical codes ( e.g. recfast 
bv ISeager et al.' (I l999l ) and IWong et all ll2008ri. RICO b y 
iFendt et al. (2009 ), RecSparse by ICrin fc Hiratal 
HyRec by |A li-Haimox id fc Hiratal (|2010al l. CosmoRec by 
IChluba fc Thomaa {,2oW ). 



2 COSMOLOGICAL MODEL 

The standard cosmological model is used. Hubble constant 
_ff as a function of redshift z is given by the following equa- 
tion: 

H{z) = HoV^A + nK{l + zY + f7„(l zY + a-e;(l + ^Y(X) 

where He, is the value of Hubble constant in present epoch, 
f^A is the vacuum-like energy density, Q.k is the energy den- 
sity related to the curvature of the Universe, Q.^ is the en- 
ergy density of non-relativistic matter, which includes con- 
tributions from CDM and baryonic matter f2m = ^cdm + 
Jlfc. firei is the energy density of relativistic matter, which 
includes contributions from photons (mainly CMBR) and 
neutrinos firei — + Qv. 

The photon energy density is related to the temperature 
of CMBR: 



(2) 



where an is the radiation constant. To is the CMBR tem- 
perature at present epoch, pc = SHq / (SttG) is the critical 
density of the Universe, c is the speed of light, G is the 
gravitational constant. 

At the present epoch the relativistic neutrino energy 
density is related with the photon energy density by the 
following formula: 



n) 



4/3 



(3) 



where Ni, is the effective number of neutrino types. 

In this paper (and current version of code) we consider 
Sltot = i^A + + + ^rei = 1, SO fiif is fixed by relation 

Qk = 1 — i^A + + f^rei)- 



The total concentrations of hydrogen and helium de- 
pend on redshift by the following formulas: 



Nh = Nho (1 + zf , Nne = Nneo (1 + zf 



(4) 



where Nho and NhbO are the values of concentrations at 
present epoch. The total concentration of the primordial 
hydrogen atoms and ions at present epoch is given by the 
following relation: 



Nho ~ — —flbNp 
niH 



(5) 



where rriH is the hydrogen atom mass, Xp is the primordial 
hydrogen mass fraction (i.e. hydrogen mass fraction after 
the primordial nucleosynthesis). 

The total concentration of the primordial helium atoms 
and ions at present epoch is given by the following relation: 



Nneo 



Pc 

rriHe 



(6) 



where mne is the helium atom mass, Yp — {1 — Xp) is the 
primordial helium mass fraction. Fractions of other elements 
are considered negligible. 

The temperature T of equilibrium radiation background 
depends on redshift according to: 



T = To{l + z) 



3 HYDROGEN RECOMBINATION 



(7) 



3.1 Main equations 

The time-dependent behaviour of hydrogen ionization frac- 
tion in the isotropic homogeneous expandin g Universe is de- 
scribe d by the follo wing kinetic equation l|Zeldovich et al.l 
(il968h : IPeebie3 ifigeS Xl: 



xhii = -CHi[aHii {Tm) N^xhii - Phi [T) exp ( - 



hVa 



xhi 



where xhii = Nhii / {Nh + Nhb), Chi is the factor by 
which the ordinary recombination rate is inhibited by the 
presence of HI Lya resonance-line radiation, anii is the to- 
tal HII— >HI recombination coefficient to the excited states 
of HI, Ne = XsNh is the free electron concentration {xe 
is the free electron fraction in common notation, see e.g. 
recfast). Phi is the total HI— ^HII ionization coefficient 
from the excited states of HI, is the Lya transition fre- 
quency, Tm is the kinetic temperature of the electron gas, 
Xhi = Nhi / {Nh + Nhc) is the neutral hydrogen fraction. 
Note that ionization fractions of ionic components are de- 
fined relative to the total number of hydrogen and helium 
atoms and ions {Nh + Nhb) while free electron fraction x^ 
is normalized to the concentration of hydrogen atoms and 
ions, Nh, as it is accepted commonly (so it is necessary for 
possible use of atlant results by other numerical codes). 

The specific form of inhibition coefficient depends on 
fine effects which are taken into account. In the current ver- 
sion (1.0) of the code only radiative feedbacks for resonant 
transitions are taken into account from the whole list of 
known (considered until now) fine effects. Other fine effects 
are planned to be included in future works. So, the inhibition 
coefficient is given by the formula: 



Chi ~ 



A'^jf + A2sls 

Phi + Alfj + A2sis 



(9) 



(8) 



where A]^jj is the total effective coefficient of np-f^ls tran- 
sitions, A2ais ~ 8.22458 s~^ is the c oefficient of 2s— » ls two- 
photon spontaneous transition (e.g. iGoldmanI (|l989l )). 

The recombination coefficient anii is given by the fol- 
lowing approximation: 



aHii {T) = Fh- 



an 



(10) 

where Fh is the hydrogen fudge factor by ISeager et al.l 
l|l999l ). Ta = T[K]/W^, and a = 4.309 • IQ-^cm^s-S 
b = —0. 6166, c = 0.6703 , d = 0.5300 are the parameters 
fitted bv lPeauignot et all (|l99ll '). 

The ionization coefficient j3Hi can be found by using 
principle of detailed balance: 



I3hi (T) = anil (T) (T) exp 



(11) 



where Qe (T) — /h^ is the partition function 

of free electrons, i^c2 is the HI c— >2 transition frequency (here 
symbol "c" denotes continuum state). 

Note that recombination and ionization coefficients in- 
cluded in the kinetic equation Q should be calculated 
at diffe rent temperatures (see e.g. lAli-Hai'moud fc Hiratal 
l|2010ah V Since the kinetics of recombination process is 
determined by the free electron distribution function, the 
recombination coefficient depends on kinetic temperature 
of free electrons Tm- The kinetics of ionization is deter- 
mined by the photon distribution function, therefore ion- 
ization coefficient depends on temperature of photons T. 
Thus the detailed balance does not take place in considered 
case, and relation has mathematical sense only and al- 
lows us to avoid direct calculation of ionization coefficient, 
(3hi (T), from integral of collisions for photons and hydro- 
gen atoms. Also exponential term in ionization part of ([8]), 
exp {—hua/ [fcsT]), should be calculated at the temperature 
of radiation. The distinction of temperatures used for cal- 
culation of ionization coefficient s in the pr esent code and in 
recfast (Seagc r ct al. (1999); Wong et al.l ( 2008 ), there the 
temperature of matter is used for this aim) leads to a little 
but important difference in the free electron fraction (see 
Fig.[3. 

The equation ((8| is solved numerically together with 
other main kinetic equations (|19p . p5|) by using second- 
order method of integration of ODE system. The relative 
deviation between results obtained by atlant and recfast 
for the period of hydrogen cosmological recombination is 
presented in Fig. [T] 



3.2 Radiative Feedbacks 

Due to a great number of fine effects the radiative feed- 
backs for resonant transitions have been chosen for includ- 
ing in the first published version of recombination code. It 
is because physics of this effect is clear (it is difficult to 
state this is true about many other fine effects) and in- 
dependently obtained resu l ts of calculation of this effect 
l|Chluba fc SunvaevI (|2010al ): [Kholupenko et all (|201Q^ ') con- 
firm each other. Inclusion of feedbacks into the cod e is based 
on formulas suggested bv lKholupenko et al.l l|2010l ). Accord- 
ing this the total effective coefficient of np-f-^ls transitions 
is: 
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Figure 1. Top panel: the relative difference Axe/xe between 
free electron fractions calculated by atlant and recfast for the 
period of hydrogen recombination. Bottom panel: (i) solid curve 
corresponds to the relative difference Axe/xe between free elec- 
tron fractions calculated by atlant and recfast [i.e. the same 
as in the top panel but in logarithmic scale], (ii) dashed curve 
corresponds to the relative difference Axe/xe between free elec- 
tron fractions calculated by atlant and recfast with corrected 
ionization rate [i.e. ionization coefficient /3[jf (T) and exponen- 
tial term exp(— hi/a/ [^sT']) in Eq. ||8ll of recfast are calculated 
at the temperature of radiation T] , (iii) dotted curve corresponds 
to the relative difference Axe/xe between free electron fractions 
calculated by atlant with modified ionization rate [i.e. ionization 
coefficient liHl(Tm) and exponential term e-xp {—hva/ [kgTm]) 
m Eq. (Ell of atlant arc calculated at the temperature of matter 
Tjn] and recfast. Dashed and dotted curves are partially over- 
lapped. 
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where 5a is given by the following expression 

En — E2 



(12) 



5a = C„ exp 



ksT 



where u„ is the frequency of 71 — > 
E,i — hfn, and coefficients Cn are: 



Cn 



rH{z'n)/rH{z)) 



(13) 

1 transition {1/2 = i^a), 
(14) 



where in turn Th is the relative overheating of Lya radiation 
(occupation number rja) in comparison with its equilibrium 
value (occupation number rj^l): 



Th = {rja/ril - l) 
and 

z'„ — ((1 + z)v„+i/v„ — 1) 



(15) 



(16) 



The set of formulas (|12)I16|) allows us to calculate ioniza- 
tion history taking into account the feedbacks for hydrogen 
71 — !> 1 resonant transitions with principal quantum num- 
bers n ^ Umax- For taking into account the feedback effect 
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Figure 2. Dashed line corresponds to the relative difference 
between atlant with fudge function 5x and base result [i.e. 
Axe/xe = &x\- Parameters of fudge function are the following: 



2.166 ■ IQ- 



1019, Azp = 180. For comparison: solid 



line corresponds to the relative difference between feedback result 
for rimax = 10 [i.e. feedbacks 11 =►...=> 2 are taken into account] 
and the base result [i.e. without feedbacks]. 



we use simple perturbation theory: at the first stage (unper- 
turbed equations) atlant calculates the relative overheating 
Vh and stores this, at the second stage atlant solves equa- 
tions perturbed by the feedback effect. Relative deviation of 
free electron fraction between perturbed (rimax ~ 10) and 
unperturbed calculations is presented in Fig. [21 



3.3 Fudge factors 

Developed recombination code allows users to use fudge fac- 
tors as well as add physical effects. This opportunity has 
been included because of the following reasons: 1) list of fine 
effects leading to 0.01 - 0.1 % corrections in the ionization 
fraction may be incomplete, and it is difficult to say how long 
this list will be in final form and when this will be achieved; 
2) considerations of some fine effects give contradictory re- 
sults due to not completely clear physical pi cture (e.g. two- 



AVong fc ScottI 



200g) ] reco il 



phot on transitions from high excited states 

rioOS'). Hir ata (2008), and Labzow skv et al 

[GrachCT&rDubrovichl |200^ and lHirata fc Forbes! li2009ll ] 
and others). 

First fudge factor is the common (see recfast) hydro- 
gen fudge factor Fh introduced in the expression (|10l) for 
hydrogen recombination coefficient anii- 

Second fudge factor is the perturbation function 5x 
modifying free electron fraction: 



(17) 



where x^*^" is the free electron fraction (normalized by the 
total concentration of hydrogen atoms and ions) being the 
final result of code running (i.e. it is value shown in the re- 



sulting file), x%'^^'^ is the solution of ODE's system describing 
ionization fractions and free electron fraction. 

Note that function 5x affects final result Xe but not 
ODE's system. This allows user to control changes of free 
electron fraction strictly by including fudge factors. 

Analyzing previous works devo ted to the fine effects 
of cosmological recombina t ion (e .g. iGrachev fc DubrovichI 
(20081): [ciiluba fc Sunvae^ l|2009b')l one may note that typ- 
ical form of corrections to the ionization fraction can be de- 
scribed by a bell-shaped function (e.g. Lorentzian or Gaus- 
sian profile or others). In this work the Lorentzian function 
has been chosen to describe uncertain deviations (until now) 
of free electron fraction from well known ODE's solution: 



1 + [(z - Zp) jl^ZpX 



(18) 



where Ay is the relative amplitude of perturbation of ion- 
ization fraction, is the redshift of perturbation maxi- 
mum, A2:p is the half- width of perturbation function at half- 
altitude. 

Result of use of fudge function is shown in the Fig. [2] 
where we have plotted t^x^jx^ = 5x as function of redshift 
z. Here we show how fudge function can mimic real fine 
corrections by means of example of feedback correction for 

^max — 10. 



4 HEII->HEI HELIUM RECOMBINATION 

The time-dependent behaviour of Hell fraction in the 
isotropic homogeneous expandin g Universe is descr i bed b y 
the following kine tic equation l|Kholupenko et al] |20o3); 
I Wong et al] j2008D ): 

XHell = ~Cpar[aparNeXHe.II — — /3parexp ( — , ) XHe.l] 

93 \ ksT I 



9g \ ^B-L 



E„ 



XHel] (19) 



where Xueii ~ NhbII / {Nh + Nhb) is the fraction of Hell 
ions relative to the total number of hydrogen and helium 
atoms and ions, Cpar is the factor by which the ordi- 
nary recombination rate is inhibited by the presence of 
Hel 2^p l^s resonance- line radiation, apar is the total 
Hell— !>HeI recombination coefficient to the excited para- 
states of Hel, subscript a denotes 2^s state of Hel atom, 
Ne is the free electron concentration, ga ~ 1 is the statis- 
tical weight of 2^s state of Hel, = 1 is the statistical 
weight of I'^s state of Hel, /Spar is the total Hel^Hell ion- 
ization coefficient from the excited para-states of Hel, Eag is 
the 2'^s — >■ l^s transition energy, Cor is the factor by which 
the ordinary recombination rate is inhibited by the presence 
of Hel 2^p l^s resonance-line radiation, Oor is the total 
Hell— >-HeI recombination coefficient to the excited ortho- 
states of Hel, subscript a' denotes state 2^s of Hel atom 
ga> = 3 is the statistical weight of 2^s state of Hel, /3or is 
the total Hel— ^Hell ionization coefficient from the excited 
ortho-states of Hel, Ea'g is the 2^s l^s transition energy, 
XHei ~ NheI / {Nh + Nhc) is the neutral helium fraction. 

The para- and ortho- recombination coefficients 
are given by widely u sed approximation formulas (e.g. 
IVerner fc Ferlandl (119961)) parame ters of which are based on 
data bv lHummer fc StorevI l|l998h ): 
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as (T) = 9s ( |r 



-1/2 



-(i-ps) 



-(i+ps) 



(20) Table 1. Parameters of the approximation of 



where subscript S t akes values 'par' or ' or', and qs, ps,Ti = 
■j^qS.iw k, T2 = 3 K (|Seager et al.1 (|l999l )') are the parameters 
of approximation. For the recombination via p ara-states we 
have Qpar = 10~^°'^*''cm^s~^, ppar = 0.711 (|Seager et al.l 
For the recombination via ortho-states we have 
= 10-^"-^°''cm^s-\ por = 0.761 (jWong et al.l ^200^ )). 
The para- and ortho- recombination and ionization co- 
efficients are related by the following formula 

/3s = ^as3.(r)expf--^V (21) 

gs \ kbT J 

where subscript "c" denotes continuum state of Hel atom, 
gc = i is the statistical weight of continuum state of 
(He'''+e~), Epar is the c — > 2^s transition energy, Eor is 
the c — >■ 2^s transition energy. 

The inhibition factor Cpar is given by the following ex- 
pression: 



Cpar 



iOb/ga) AbgPtg eXp (-Eta/kBT) + Agg 

Ppar + (gb/ga) AbgPbg exp {~Eba/kBT) + Aag 



(22) 



where Atg is the Einstein coefBcient [s~^] of 2^p o l^s spon- 
taneous transitions, Ptg is the probability of the uncompen- 
sated 2^p l^s transitions. Eta is the 2^p — > 2^s transition 
energy, Aag is the coefficient of two-photon 2^s — >■ l^s spon- 
taneous decay. 

The inhibition factor Cor is given by the following ex- 
pression 



Cor 



{gv/ga') AygPyg exp {-Eyg, /kBT) 
Por + {gb' /ga') AygPi,,g exp {- Ef a' / k bT) 



(23) 



where subscript b' denotes state 2'^p of Hel, gy = 9 is the 
statistical weight of 2^p state of Hel, where Abg is the Ein- 
stein coefficient [s~^] of 2^p -f-)- l^s spontaneous transitions, 
Ph'g is the probability of the uncompensated 2^p — > l^s 
transitions, E^'a' is the 2^p — 2'^s transition energy. 

The probabilities Ptg and P^'g take into account the 
escape of Hel resonant photons from the line profiles due to 
the cosmological expansion and destruction of these photons 
by neutral hydrogen. They can be found by the following 
formula: 

Pf.-P}, + P"9 (24) 
where PJg is approximately given by the following: 

^(l-exp(-rHe,/)) (25) 



-1 

THe. 



where 7 is the ratio of the helium and hydrogen absorption 
coefficients at the central frequency of the /—>■(/ line (here 
symbol f — b or b' depending on what transition is con- 
sidered), THe.f is the Sobolev optical depth. The value 7 is 
given by the following relation: 



7 = 



igf/gs)AfgNHelC^ 



(26) 



where an is the ionization cross-section of hydrogen ground 
state, parameter /^VD,f ~ ^fg\/'^kBT / (niHeC^) is the 
Doppler line width. 

The optical depth THe.f is: 

rHej = gfAfgNneic''/ {ggSnHujg) (27) 



Range of 7 




P 


q 


s; 7 ^ 5 ■ 10^ 




0.66 


0.9 


5 • 10^ < 7 s; 5 ■ 


10* 


0.515 


0.94 


5 • 10"' < 7 s£ 5 ■ 


10^ 


0.416 


0.96 


5 • 10^ < 7 




0.36 


0.97 



The value Pfg is given by the following: 

Pfg ^pS + Pg+Pr (28) 
where in turn P§ approximately is: 

pS = {1 + pj")-' (29) 

where parameters p, q are given in the Tab. [T] 

The value Pq is given by the following approximate 
formula: 



pH 



(8Aa)'' -TT --7 
[1 -f exp (-1.07 In (s -I- 1.5) - 0.45)] (s + 1.28)"^''^ 



(30) 



where a — Vf/AnAvoj is the Voigt parameter (here F/ 
is the natural line width), A is the single-scattering albedo 
in / ^ g line, and par ameter s is given by the following 
relation icrachevl (119881 )): 

s = 2-^/2^-^/^(1 - \)\-'/\^'^^'/^ - 1/4 (31) 
The single-scattering albedo A can be calculated by: 

where Rf = "^^Rf^n is the total coefficient of radia- 
tive transitions (including induced transitions, while colli- 
sion transitions are considered negligible) from / state to 
excited states of Hel atom. 

The value P^ is approximately: 



2(1-A) 



arctan (/(a. A, 7)) ) (33) 



(1 - A)77r3/2 

where /(a. A, 7) is defined by the expression 

/(a. A, 7) ~ 0.69 • (1 - A)-'/2a'/* ("7)^''"' (34) 

The equation p9|l is solved numerically together with other 
main kinetic equations ([8]) and (|35|l to determine free elec- 
tron fraction Xe- The relative deviation of free electron frac- 
tion between the results calculated by using atlant and rec- 
fast for the period of HeH— >-HeI recombination is presented 
in Fig. El 

5 HEIII->HEII HELIUM RECOMBINATION 

In the difference with recfast we use non-equilibrium (i.e. 
kinetic) approach for consideration of Helll— >-HeII helium 
recombination. The time-dependent behaviour HeHI frac- 
tion in the isotropic homogeneous expanding Universe is de- 
scribed by the following kinetic equation: 



XHelll = —CHeIl[0!HeIII NeXHelll 
hUHell, 21 



PhbII exp 



fcsT 



XHell 



(35) 




' 1600 1800 2000 2200 2400 2600 2800 3000 
z 

Figure 3. The relative difference Axe/x^ between free elec- 
tron fractions calculated by atlant and recfast for the period 
of Hell— >-HeI recombination. 




z 



Figure 4. The relative difference Axe/xe between free elec- 
tron fractions calculated by atlant and recfast for the period 
of HcIII— >HeII recombination. 



where XHeiii = Nh^iii/ [Nh + Nhb), ChsII is the factor 
by which the ordinary recombination rate is inhibited by the 
presence of Hell 2p — >■ Is resonance-hne radiation, aneiii 
is the total Helll— >HeII recombination coefficient to the ex- 
cited states of Hell, PhbII is the total Helf^HeHI ioniza- 
tion coefficient from the excited states of Heff, UHe.n,2i is 
the Heff 2p — > la transition frequency. 

The inhibition coefficient is given by the expression: 



ChbII = 



A 



HeII.2pls 



I3l 



A 



HeII.2sls 



(36) 



where A^u^ji 2pis is the effective coefficient of Heff 2p^fs 
transitions, AHeii,2si3 is the coefficient of Heff 2s— >-fs two- 
photon spontaneous transition. 

Since Hefl is the hydrogenic ion the recombination coef- 
ficient a Heiii can be found by us ing simple scaling relation 
(see e.g. IVerner fc Ferlandl (|l996l )): 



az (T) = Zoii {T/Z'^ 



(37) 



where Z is the nuclear charge for the hydrogenic ions. In 
considered case Eq. (p7)) yields aneiii (T) = 2aHii [T /A), 
where anii is taken from (|f 0|) without hydrogen fudge factor 
Fh- 

The ionization coefficient PhbII can be found by using 
principle of detailed balance: 



Pnell (T) = OlHelll (T) (T) CXp 



(38) 



i'Heii,c2 is the Heff transition frequency. 

Note that in Eq (|35|l the recombination and ionization 
coefficients are calculating at the same temperature T. This 
is because the matter temperature Tm is very close to the 
radiation one T (relative deviation is less than fO~^) dur- 
ing HeHI— >^Helf recombination and calculation at different 
temperatures has no sense at required level of accuracy. 

The effective coefficient A^i]^ij 2pis of Hell 2pof s tran- 



sitions due to escape of Hell 2p-f-^fs resonant photons from 
the line profile because of cosmological redshift is given by 
the following formula: 



A 



HeII,2pls 



(39) 



The value AHeii.2si3 is found from charge scaling for the 
hydroge nic ions Az.2sis ~ Z^ A2s \b (see e.g. , Shapiro fc Breit 
(195^, IZon fc RapoportI !\196^ . iNussbaumer fc Schmutl 
(298J)). In the considered case this gives us AHeii,2sis ~ 



The equation psp is solved numerically together with 
other main kinetic equations and (|19p to determine free 
electron fraction Xe- The relative deviation between results 
by atlant and recfast for the period of Helll— >-HeIf recom- 
bination is presented in Fig. [J] 



6 EVOLUTION OF MATTER TEMPERATURE 

The behaviour of the matter temperature in the isotropic ho- 
mogeneous ex panding Unive rs e is described by the following 
equation (e.g. ]Peeblej l| 19681 ) : IScott fc Mosd i|2009^ ): 



Tfji 



SaraRT" 



Sm^c 1 + Xe + Nhbo/Nho 



{T - Tm) - 2HTm{4:0) 



where ar is the Thomson scattering cross section. 

Defining relative deviation, 5t, of the matter temper- 
ature from the radiation temperature via Tm = T (1 — St) 
and substituting this into (|4Up one can obtain: 

5t = - {Rt + H)5t + H (41) 

where Rt is the rate of energy transfer between matter and 
radiation via Compton scattering: 

SaraRT* Xe 



Rt 



3meC l + Xe+NHeo/NE 



(42) 




z 



Figure 5. The corrections, Sj-, to tiie matter temperature, Tm, 
relative to the radiation temperature, T, as functions of redshift z: 
dashed curve corresponds to the correction, (5y g i dashed-dotted 
curve corresponds to |5t,i1 {5t,i < 0), and solid curve corre- 
sponds to ((5t,0 + <^T,l)- 

To solve equation (|4ip we applied a perturbation ap- 
proach. In the early stages of the Universe history (z > 200) 
the rate of energy transfer between matter and radiation via 
Compton scattering is much larger than the rate of the tem- 
perature change (the latter is about Hubble expansion rate 
H), i.e. St / (RtSt) ^ 1. Thus we can use expansion of the 
solution over this smallness: 

oo 

St = Y, ^T,i (43) 

where zeroth-order approximation is deter mined a s qua - 
sistacionary solution o f Eg. (|4T|| (see also iHiratal l|2008l ): 
lAli-Haimoud fc Hiratal (|2010bfr i7 

5t,o = {Rt/H + 1)-^ (44) 

and the next members of expansion (|43l) are related by the 
following equation: 

ST,r+l = - {Rt + H)-^ 5T,i (45) 

In the present version of the code we keep only two first 
corrections St,o and St,i in Eq. (|43p . 

For the period z < 200 the quasistationary condition 
is violated, so expansion H43p loses convergence and cannot 
represent the solution of Eq. H4ip . In this redshift range we 
use the following dependence of matter temperature on red- 
shift: 

T,n = r^^= f ^^^V (46) 

where T!^'^ is the matter temperature from (|43p at z^ec (in- 
dex "dec" means "decoupling") which is considered as the 
moment of decoupling of the matter temperature from radia- 
tion the one (at value z^ec we should make join of solutions). 
In present version of code we determine Zdec from condition 



^ \ ' \ ' r 




z 



Figure 6. The relative difference Axe/xe of the free electron 
fractions as a function of redshift z: dotted curve corresponds to 
the difference, Axe/xe, between results for (Sy = (i.e. Tm = T) 
by atlant and recfast corrected (for the details of corrections, see 
section FS.ll of the present paper), dashed curve corresponds to the 
difference, Axe/xe, between results for St = <5t,o by atlant and 
corrected recfast, and solid curve corresponds to the difference, 
Axe/xe, between results for St = {pT,o + ^T,i) by atlant and 
corrected recfast. 

5t,o = 0.5. Equation (|46|l corresponds to the equation of 
state of non-relativistic matter. 

Results of calculations of 5t,o and 5t,\ are presented in 
Fig. [H The influence of taking into account different approx- 
imations for 5t [depending on the number of kept members 
of expansion (|43|) ] on the free electron fraction is shown in 
Fig. El 



7 VARIATION OF THE FUNDAMENTAL 
CONSTANTS 

Today the opportunity to vary of the fundamental constants 
becomes essential at the analysis of CMBR anisotropy (e.g. 
IScoccola et ai] l|2008l )) and we decided to include this in our 
code. Thus the current version of the code allows user to see 
how recombination occurs at different values of the funda- 
mental constants (this means that changes of fundamental 
constants lead to corresponding changes of derived physi- 
cal values, e.g. ionization energies of atoms, Thomson cross 
section and others). 

Since some of used physical values are given only nu- 
merically (e.g. level energies and transition probabilities for 
Hel atom, see www.nist.gov) we use simple scalings to take 
into account the influence of variation of fundamental con- 
stants on these values. These scalings are the following: 

1) For the level energies: 

E = eJ^X(^)(±Y (47) 
\est J yrrie.st ) \nst ) 

2) For the one-photon transition coefficients: 




1000 



2000 



3000 



Figure 7. Free electron fractions Xe as functions of redsiiift z 
at different values of the fine-structure constant afi„e- dasiied 
curve corresponds to afi^^/a^^^^^ = 0.95, solid curve corre- 
sponds to Oifine/ci'fing = 1> daslied-dotted curve corresponds to 



A'' = Al 



, 10 , 

e \ / rrie 



h 



^Cst J \me,st J \h, 

3) For the two-photon transition coefficients 



A^^ = AZ7 — 



, 16 , 

e \ / me 



^est J \me.st 
4) For the recombination coefficients 



a — ast 



a / \ -3/2 / t \ -1 / \ -3 



est 



me, St 



hst 



Cst 



(48) 



(49) 



(50) 



In the expressions (|47|) - (|50| l subscript "st" denotes current 
values of physical quantities while symbols without subscript 
denote values which user suggests to be valid during the 
recombination epoch. 

The influence of variations of the fine structure constant 
a fine (via varying elementary charge) on the free electron 
fraction is shown in Fig. [T] 



8 RESULTS 

The main numerical results of this paper are the relative 
differences of free electron fractions Axe/xe calculated by 
atlant and recfast (figs [T] [31 13| . In the top panel of Fig. [T] 
the relative difference Axe/xe between result of atlant an d 
current version of recfast (version 1.5. IWoim et al.l (|2008l )) 
for the hydrogen recombination epoch. The first significant 
difference appears in the period z = 1550 — 1650. This is 
because at early epochs z > 1550 recfast considers re- 
combination as quasi-equilibrium process (according to Saha 
formula) and sets initial condition for the kinetic equation 
correspondingly, while atlant considers recombination as a 



non-equilibrium (kinetic) process in the whole range of red- 
shifts and sets an initial conditions corresponding to the 
fully ionized plasma at the moment given by user (in our 
case Zbegin ~ 8000). The sharp peak of Axe/xe in the pe- 
riod 2 — 1550 — 1650 (majdmal value is about 5 ■ lO^**) 
is due to transition from recombination according to Saha 
formula to non-equilibrium recombination which occurs in 
recfast at free electron fraction Xe = 0.985 (that corre- 
sponds to z ~ 1555 in considered case). There is a simple 
method of modelling of CMBR s pectral distortion due t o 
cosmological recombination (e.g. [ Bern shtein et al.l (| 19771 ): 



iDubrovich fc Stolvarovl (| 19951 . [igg^HBurEin (200i)). This 
method is based on the formalism of matrix of efficiency 
of radiative transitions (ERT-matrix) and it does not de- 
mand direct calculation of atomic level populations. On the 
other hand this method demands use of the following val- 
ues dxHii/dz, dxHeii / dz, and dxHeiii/dz to determine the 
rates of resonance photon emission and correspondingly the 
shape of cosmological recombination lines. Thus for using 
ERT-matrix method the sharp numerical features may be 
significant defects of cosmological recombination modelling. 
The second difference of results by atlant and recfast ap- 
pears in the period z = 300 — 1000. It arises due to different 
structure of ionization items of hydrogen kinetic equation 
((8)l in atlant and recfast: in recfast the ionization coeffi- 
cient Phi and Boltzmann exponential term are calculated at 
the temperature of matter while in atlant at the radiation 
temperature. Most part of this difference is due to distinc- 
tion of ionization coefficients including in the denominator of 
inhibition coefficient Chi- The maximum of Axe/xe in men- 
tioned range of redshifts is about 0.27% at z ~ 770. This is 
little but maybe important difference in the context of fu- 
ture analysis of Planck data. The third difference appears at 
low redshifts [z < 300). It is due to the different approaches 
to evaluating matter temperature in recfast (where ODE 
for temperature is solved) and in atlant (where perturba- 
tion theory for temperature estimate is used) . In the bottom 
panel of Fig. [T] we present three curves. Solid curve is the 
same as in the top panel but in a logarithmic scale. Dashed 
curve corresponds to the relative difference Axe/xe between 
free electron fractions calculated by atlant and recfast with 
corrected ionization rate [i.e. ionization coefficient Phi (T) 
and exponential term exp {—hi/ci/ [ksT]) in Eq. ((8)l of the 
recfast model are calculated at the temperature of radi- 
ation, T]. Dotted curve corresponds to the relative differ- 
ence Axe/xe between free electron fractions calculated by 
atlant with modified ionization rate [i.e. ionization coeffi- 
cient Phi (Tm) and exponential term exp { — hvc,/ [ksTm]) in 
Eq. ([8]) of the atlant model are calculated at the temper- 
ature of matter Tm] and recfast. These curves show that 
in the frame of identical physical models the accordance of 
results by atlant and recfast is wholly satisfactory. Resid- 
ual difference which does not exceed 5 ■ 10~^ for the period 
z = 400 — 1500 can be explained by a little different values of 
physical constants and different integration methods which 
have been used in atlant and recfast. 

In Fig. [5] we demonstrate how feedbacks for resonant 
transitions of hydrogen affect free electron fraction (solid 
curve) and show an example of using fudge function, 5x 
(dashed curve). The change due to feedbacks has the maxi- 
mum atout_0;2166%_atredsM z ~ 1019 in full accordance 
with lChluba fc SunvaevI l|2010ah . This calculation has been 
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done for Umax ~ 10. Parameters of the fudge function have 
been choosen as foUowing: Ap = 2.166- 10"'^, Zp — 1019, and 
Azp = 180. This is to show how use of the fudge function 
allows us to imitate the influence of real physical effects on 
the free electron fraction. 

In Fig.[3]the relative difference of the results by atlant 
and recfast for the epoch of Hell— >-HeI recombination is 
presented. The breaks in the range z = 2700 — 2900 and 
break at z ~ 2300 are the numerical features connected 
with switch in recfast code. The maximum of difference 
is about -0.23% at z ~ 1950. The difference between re- 
sults by atlant and recfast for the epoch of Hell— >HeI 
recombination arises because different approximations for 
estimate of effective eacape probability are used: in at- 
lant the approach based on analytical consideration of res- 
onance radiation trans fer in the pres ence of continuum ab- 
sorption is used (Kh olupenko et al.l (^2008)) while recfast 
uses the si mple approx i matio n with an additional fudge 
factor bHe (|Wong et al.1 (|2008l lV It leads to the numerical 
differences because used approximations are not identical. 
Also there is physical rea son for divergence of these ap- 
proximations: IWong et al.l (|2008 ) have determined bne from 
the best agreement with result bv lSwitzer fc Hiratal (|2008al ') 
who took into account not only hydrogen continuum opac- 
ity but also feedbacks for resonant transitions of Hel and 
foun d this effect is about 0.46% at z ~ 2045. Later works 
fe.g. lChluba fc SunvaevI (|2010al )') show that this feedback ef- 
fect is about 0.17% at z ~ 2300, i.e. negligible for the most 
of modern problems connected with cosmological recombi- 
nation (e.g. CMBR anisotropy analysis). Thus we do not 
take helium feedbacks into account in the present version 
of the code. One more important point of calculations of 
effective escape probability for Hel is the estimate of num- 
ber of neutral hydrogen at the high redshifts z > 1600 
(i.e. when hydrogen ionization fraction close to unity [e.g. 
greater than 0.985]): recfast uses Saha formula for this aim, 
while in atlant the kinetic equation is solved in the whole 
range of redshifts to determine the fractions of all considered 
plasma components. This explains the discrepancy between 
results of the present pa per and Kholupenko et al. (200i) 
where approximation by IWong et al.l (|2008l ) has also been 
investigated by using earlier version of our code (i.e. Wong- 
Moss-Scott approximation has been calculated, but number 
of neutral hydrogen has been calculated from kinetic equa- 
tion). 

In Fig. 0] the relative difference of the results by atlant 
and recfast for the epoch of Helll— >-HeII recombination is 
presented. The break at z = 5000 is the numerical feature 
connected with switches in recfast code again. In the range 
of redshifts z — 5000 — 7000 the difference between the re- 
sults by atlant and recfast arises because Helll^Hell re- 
combination is treated by atlant in non-equilibrium way 
while recfast treats this according to Saha formula. The 
maximum of difference is about 0.225% at z ~ 5790. Such 
difference is too small to be a reason of any observational 
effects at the current level of experimental accuracy, but 
non-equilibrium treatment of Helll— >HeII allows us to avoid 
artificial numerical features (that is important for the mod- 
elling of CMBR spectral distortion arising from Helll— >-HeII 
recombination) . 

In Fig. [5] the deviation of the matter temperature 
from the radiation temperature is presented. The dashed 



curve corresponds to the correction St.o, dashed-dotted 
curve corresponds to |(5t,i1i a-nd solid curve corresponds to 
{St,o + 5t,i)- One can see that difference between the ra- 
diation and matter temperature is less than 10~* down to 
z ~ 940 and less than 10~^ down to z ~ 790. Taking into 
account that Tm is included in kinetic equation only 
via recombination coefficient ohii (Tm) one can expect the 
similar changes in hydrogen ionization fraction at transition 
from approximation Tm = T to more exact estimate of Tm- 
From accuracy point of view the correction 5t,i becomes 
important beginning from z ~ 480 where/when it achieves 
values about 10""^. At the low z {z < 300) the correction St,i 
is larger than 20% of St.o and should be taken into account 
for correct determination of residual ionization fraction as 
possible. In Fig. [S] we show how including of ST,i's affects the 
free electron fraction. We have compared our results with the 
results by recfast corrected in the following way: ionization 
items of Eq. ((S)), i.e. ionization coefficient /3hi and exponen- 
tial term exp { — hva/ [fesTm]) are calculated at temperature 
of radiation T. The dotted curve corresponds to the relative 
difference Ax^jx^ of the result by atlant for approximation 
Tm ~ T (i.e. 5t = 0) and the mentioned result by recfast, 
the dashed one corresponds to Ax^/x^ at 5t ~ St.o, and 
the solid one does Axe/xe at 5t ~ {St,o + St,i)- From Fig. 
[6] one can see that rough approximation Tm = T is quite 
valid (at required level of accuracy) down to z ~ 670 where 
Axe/xe achieves values about 10^"^. Thus using this rough 
approximation biases the results not so dramatically as it 
seems when this approximation is used in recfast. More ad- 
vanced approximations are valid [Axe/xe ^ 10"'^) down to 
z ~ 310 {St = St,o) and z ~ 250 {St = {St,o + St,i)). 

Fig. [7] illustrates one of the additional opportunities of 
atlant: here the changes of free electron fraction at varia- 
tion of the fine-structure constant are presented. Presented 
results are obtained for the following values of the fine- 
structure constant afi„e'- O.QSa'ji^^, a^^^^, and 1.05a°i„g 
(where a/i„e is the current value of the fine-structure con- 
stant). Obtain e d res ults are very similar to the results by 
IScoccola et al.1 (|2008l ). 
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